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Abstract 

The exponential of an A x A matrix can always be expressed as a matrix polynomial of order A — 1. In 
particular, a general group element for the fundamental representation of SU{N) can be expressed as a matrix 
polynomial of order A — 1 in a traceless A x A hermitian generating matrix, with polynomial coefficients consisting 
of elementary trigonometric fnnctions dependent on A — 2 invariants in addition to the group parameter. These 
invariants are just angles determined by the direction of a real A-vector whose components are the eigenvalnes of 
the hermitian matrix. Eqnivalently, the eigenvalues are given by projecting the vertices of an (A — l)-simplex 
onto a particular axis passing through the center of the simplex. The orientation of the simplex relative to this 
axis determines the angular invariants and hence the real eigenvalues of the matrix. 


“Let us revisit Euclid. Let us discover for ourselves a few of the newer results. Perhaps we may be able to 
recapture some of the wonder and awe that our first contact with geometry aroused.” — H S M Coxeter 


Curtright and Zachos (CZ) wrote a brief summary [T] of essential elementary features for the triplet representation 
[2] of SU{3), thereby distilling some older results [51I11[S]. Here, I show how their main results may be extended to 
almost any N x N matrix, thereby embellishing [aiisiiT]. A matrix polynomial form for an exponentiated matrix 
can always be expressed succinctly in terms of a “response function” that encodes the eigenvalues of the matrix. 

In particular, I show how the CZ results may be extended to the fundamental representation of SU{N) for any N. 
I show that a polynomial form for any group element in the fundamental representation can be expressed in terms of 
a response function that encodes the real eigenvalues of the hermitian matrix that generates the group element. In 
addition, I provide a clear geometrical picture of the relevant group invariants in terms of the elementary properties 
of an (A — l)-simplex [8], as a simple generalization of Viete’s venerable results for the real roots of a cubic equation 
[5]. I give specific results for SU (A < 5). 

The exponential of an A x A matrix M can be written as a matrix polynomial as a consequence of the Cayley- 
Hamilton theorem m, or more directly, as the result of the Lagrange-Sylvester projection matrix method M- A 
reasonably compact polynomial form is given by 

Af-l 

exp (itM) = ^ M'^En (t) , 

n—0 

N — l — n / 7 \ N —l—n—Tfi 

Enit)= 

m—0 ^ 

The invariant functions appearing in the En coefficients are given by symmetric polynomials, 

Am [A] = ^ I (3) 

l<A:i <^2 

where for A; = 1, • • • , A are the eigenvalues of M, and by various derivatives of the response function for the 
matrix, F (t). The latter is defined in terms of the characteristic function C {z) as 

- C iz) = det {zl - M) , (4) 

where I have assumed the eigenvalues are non-degenerate. For degenerate eigenvalues, appropriate limits must be 
taken. 

Alternatively, the symmetric polynomials can be written as invariant traces, hence computed directly from M 
without knowing the individual eigenvalues |12) . 
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This form of the invariants follows from the generating function for the polynomials as given by a determinant, 

N 

det (/ + tM)=Y^ ■ (7) 

m—0 

The result © may also be derived in a straightforward way just by taking the inverse Laplace transform of the 
resolvent for a general complex N x N matrix M, as expressed by 

1 ^-1 n 

7VTF =!:«■■«.(»). A. (») = ‘ • <*> 

n=0 ^ ' 

where Trunc [/ (s)] = X;™=o (0) is the fcth-order Taylor polynomial for the function /. This result for 

k _ 

the resolvent is well-known [nnnin]. In this approach the response function is encountered as 

/ Atz 


where the counter-clockwise contour integral encloses all the eigenvalues of M. 

To be more specific I now give some results appropriate to the group SU (N), for various N, that evince most 
of the general features. In these specific cases I is always the N x N unit matrix and iL is a traceless, hermitian 
matrix, so the eigenvalues are all real. In addition, Fm (t) is the response function for the N x N matrix H, as 
defined above. The corresponding characteristic function and its derivative are 
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Again, I will assume the eigenvalues are non-degenerate. Otherwise, for the expressions to follow appropriate limits 
must be taken. 

For the SU (2) doublet, as is well-known. 


exp {itH) = 


H -il 


dt 


F2{t) 


For the SU (3) triplet, as given in [T], 

exp {itH) = 

For the SU (4) quartet. 
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And finally, for the SU (5) quintet. 
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Etc. Note how the form of the result for SU (N) can be immediately obtained from that for SU {N + 1) by discarding 
the unit matrix term for the latter, by decrementing the exponents of the remaining matrix powers of H (but do not 
change the exponents inside invariant traces), and by replacing EN+i —>■ F^. 

Conversely, the result for SU (N) can be obtained from that for SU {N — 1) just by replacing F^-i —>■ EN, by 
incrementing the powers of H in the SU {N — 1) expression, and finally by adding the appropriate unit matrix term. 
That unit matrix term for SU (TV) can always be expressed as a series of derivatives of Fn with “Viete coefficients” 
given by invariant traces, as described above, namely. 
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The sum here can also be written as an inverse Laplace transform: 
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The explicit form of the Ik invariants for all k are given in ([ 6 ]) as traced powers of H. For the specific examples 
it suffices to note: Iq = 1 , Xi = tr(iJ), I 2 = (ti {H))^ — tr(iL^), I 3 = (tr(7L))^ — 3tr (iL) tr (iL^) + 2 tr(iL^), 
and I 4 = (tr — 6 (tr (iV))^ tr (H^) + 8 tr (H) tr + 3 (tr — 6 tr {H^), where tr (H) = 0 for SU {N) 

generators. 

Thus a compact form of the fundamental exponential polynomial for any SU (N) can be obtained by using ()15p 
to construct sequentially the hierarchy of polynomials for SU {M < N), starting from the trivial result / for SU (1) 
[18) . However, more work is required to obtain explicit results for the response functions since the eigenvalues for a 
generic N x N matrix, whether hermitian or not, pose some challenges m, especially in the limit of large TV. 

Nevertheless, as an explicit example, an element of an SU (2) subgroup of SU (TV) is always manageable in closed 
form, where TV = 2j +1 fixes the spin of the embedded SU (2) representation. The standard choice for the generator 
n ■~S of a rotation about axis h gives a characteristic function 


2j 

C(A) = det(A/-h.:/) = f[(A-(j-A)) = ^E±iJ:E . (17) 

k=0 ^ 

The trace norm of this standard generator is given by a product of the quadratic su (2) Casimir, j (j + 1), and the 
matrix rank, TV = 2j + 1, namely. 


tr 


2l 1 

(h • 7 ) =^{j -kf = - j (j + 1) X (2j + 1) 
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Traces of all odd powers of h • J vanish, of course, since the eigenvalues are symmetrically distributed about zero. 
Traces of all even powers are easily computed, as given by the following generating function (i.e. the character of 
the group element, upon replacing x ^ iO where 9 is the rotation angle): 
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Note that 


2j + l 


tr 


= (2j + 1) (^1 + y j (j + 1) + ^ i (j + 1) X (3j (j + 1) - 1) + O (7)^ 

is always a kt\\ order polynomial in the quadratic Casimir. 


More importantly for the exponentiated matrix, the response function for the embedded SU (2) spin j represen¬ 
tation is 

Fn=2,+i (9) - ^ _ ^) - 7^ sm (9/2) . (20) 
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Derivatives of this response function are readily evaluated to obtain a matrix polynomial expansion for exp ^i9n ’ j 
although it takes quite a bit of additional work to reduce the resulting polynomial to the compact form in HD Hi]. 
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For generic SU (N) generators the eigenvalues and the response function cannot be so easily determined as was 
the case for the embedded SU (2) example. However, there is an elegant geometrical way to parameterize the 
eigenvalues that generalizes the well-known picture developed by Viete in the 16th century for cubic equations with 
three real roots [5]. Indeed, the discussion to follow is easily adapted to provide a visualization of the solution to 
any TVth order equation with only real roots, although it does not immediately give an explicit solution for the roots 
in terms of traces and/or determinants. 

For SU (iV) the N real eigenvalues of the traceless hermitian generator H may be envisioned as the components of a 
Euclidean vector, A G E^r. This vector of eigenvalues lies in an (TV — l)-dimensional hyperplane that passes through 
the origin of the Euclidean space, with normal T/ = (1,1, • • • , 1) G E^v, as defined by 0 = T/ • A = J2k=i = tr (H). 
The vector A is also a point on a sphere Sn-i C Ejv, with the non-vanishing radius of the sphere defined by the 
matrix invariant = A • A = ~ intersection of this hyperplane and (TV — l)-sphere define 

a particular (TV — 2)-sphere in the eigenvalue space, Sn -2 C Ejv. The eigenvalues of H therefore comprise the 
components of a vector A that is a point on this Sn- 2 , and for a given radius r the vector A may be completely 
specified by the standard N — 2 spherical polar angles parameterizing Sn- 2 - 

For example, for SU (3) the vector of eigenvalues is a point on the circle defined by the intersection of the 
aforementioned plane and 2-sphere. The components of A are specified by the radius of the circle and by a single 
angle 9, as noted long ago by Viete for cubic equations with three real roots. Thus 

3 

Afe =-\/| rcos (0-I-27rfc/3) for fc = l,2,3, where A^ = cos (30) = det (TV) . (21) 

fc=i 


Perhaps the only surprising features apparent in (ICT) are the facts that the radius of the circle has been rescaled [20] 
and that the eigenvalues may also be viewed as the projections of three points equally spaced on the circle. That 
is to say, the eigenvalues are given by projecting along a particular axis the vertices of an equilateral triangle (i.e. a 
2-simplex) circumscribed by a circle, where the axis of interest is a diameter of the circle. The angle 9 is determined 
by the orientation of the triangle with respect to this axis. 

Perhaps even more surprisingly, this geometrical picture generalizes completely to any TV. The eigenvalues of any 
generator for the fundamental representation of SU (TV) are given by projecting the TV vertices of an (TV — l)-simplex 
onto a particular axis that passes through the center of the simplex. The orientation of the simplex with respect to 
the axis determines the TV — 2 angles that are needed, in addition to the radius r, to specify the components of A . 
I now explain in detail how these geometrical facts are established, if they are not already obvious. 

The components of A are the eigenvalues of H. Each of these components is given by a Euclidean inner product 
with one of the unit vectors that define the standard orthonormal basis for Ejv. Thus 


Afc=efe-A, where ei = (1,0, • • • , 0) , 62 = (0,1, • • • , 0) , 


Cat = (0, • • • ,0,1) 


( 22 ) 


On the other hand, these TV unit vectors are the vertices of the TV-cell that defines the standard (TV — l)-simplex, 
embedded in Eat. This standard simplex lies in a hyperplane “parallel” to the one that contains the vector of 
eigenvalues, with a normal again given by the aforementioned it. Moreover, since it ■ X =0, the inner products 
with A are unchanged if each of the Ck is translated by adding any amount of the normal it. Thus we may rigidly 
translate the standard simplex so that it too is centered on the origin and lies in the hyperplane containing the 
eigenvalues. Upon doing so, we obtain a simplex with vertices Ck — , where the eigenvalues are now given by 

inner products Xk = (ck ~ ^ ■ More explicitly, 
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In addition, we may rescale these vertices to obtain a simplex that also lies on the (TV — 2)-sphere that contains the 
vector of eigenvalues. Thus take 

1 . 


= 




TV 


TV- I 

One readily checks that these TV vectors are indeed the vertices of a simplex, with 
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This simplex corresponds to the array of weight vectors m for the fundamental representation of SU (TV). The 
eigenvalues are now obtained by projecting the vertices of this final (TV — l)-simplex onto an axis defined by the 
direction of A , i.e. by projectir^ the vertices onto a particular diameter of the {N — 2)-sphere. 

Thus, with unit vector e = A /r, 

>'k = . (26) 

But note, what matters for the eigenvalues is the relative orientation of the simplex with respect to the axis. Thus 
we may produce the same eigenvalues if we choose the axis to be in any convenient direction, and appropriately 
orient an equivalent simplex with respect to the chosen axis. 

For example, for the eigenvalues of the SU (3) generator, as given in (|2T]) . we may choose the axis as the horizontal 
on the plane, with 9 the planar angle measured counterclockwise from the horizontal. We may then visualize the 
eigenvalues as projections onto the horizontal axis of the vertices of an equilateral triangle in the plane — a 2 -simplex 

— those vertices being located at angles 9 and 0 ± 27r/3, and at a distance r from the origin. 

For other examples, consider SU (4) and SU (5), and note that the freedom to re-orient jointly the axis and the 
simplex, as if they comprised a single rigid body, leads to apparently different — but nevertheless equivalent — 
parameterizations of the eigenvalues in terms of the standard spherical polar angles for Sn- 2 - 

For SU (4) the vector of eigenvalues is a point on a 2-sphere. The four eigenvalues, i.e. the components of the 
four-vector A , may be specified in terms of two polar angles. 


Ai = r 
A 2 = r 
A 3 = r 


A 4 = r 
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These are just projections of the vertices of a regular tetrahedron onto an axis through its center, whose orientation 
is left as an exercise for the reader [22], but which is easily visualized in three Euclidean dimensions. 



A tetrahedron pierced by an axis through its center. 

For SU (5) the vector of eigenvalues is a point on a 3-sphere. The five eigenvalues, i.e. the components of the 
five-vector A , may be specified in terms of three polar angles. 
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These are just projections of the vertices of a pentatope onto an axis through its center. In this case I have chosen 
things so that the relative orientations of the axis and the five vertices are easily determined in terms of the angular 
parameterization, even if the entire geometry is not so easily visualized as in the SU (4) example. 

Finally, just as the eigenvalues are readily expressed in terms of the radius and the various angles, so too are the 
latter expressed in terms of combinations of the eigenvalues. But it is not easy — there is no royal road to find the 
roots of quintic and higher order equations — to express the eigenvalues, or the corresponding angles, in terms of 
the invariant traces or the determinant of H for SU {N > 4). 

For example, for SU (4) as parameterized above, 


4 4 


tr (-ff) = X! ® ’ tr {H^) = ^ A^ = , 

k^l k^l 
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(36) 
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tr ^ sin 9 sin [29) cos (20) , 

k^l 
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det {H)^l[Xk = l ((tr (ij 2))2 _ 2 tr (H^)) 



k=l 


= (l + (2 sin^ 4’ ~ ^) sin^ d) (l + (2 cos^ — 3) sin^ 9) . (38) 

This leads to equations cubic and quadratic in sin^ 9 and sin^ ip that can be solved in closed form in terms of 
elementary functions of the traces of H^, and H'^ or det(7J). For detailed discussions on how to solve the 
resulting equations, see [61123]. 

For SU (5) as parameterized above, it is straightforward to check that 


tr {H^) = 




{cosip) 





(39) 

etc. However, in this case I can not find an elementary expression of the angles, or eigenvalues, in terms of invariant 
traces of powers of H. But then, neither can anyone else, as is well-known [24] . 

Nonetheless, at the very least I hope the constructions described here for any N may serve as elementary examples 
that illustrate the elegant interplay between group theory and geometry. As H S M Coxeter once said [5], 


“ ... in four or more dimensions, we can never fully comprehend them by direct observation. In attempting to do 
so, however, we seem to peep through a chink in the wall of our physical limitations, into a new world of dazzling 
beauty. Such an escape from the turbulence of ordinary life will perhaps help to keep us sane.” 


Acknowledgement I thank Professor Curtright for focusing my attention on [T], and for suggesting that the 
simple geometrical construction so succinctly described therein would generalize to any SU {N). 
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